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Abstract 
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1 Introduction 



In the scalar theory of gravitation by Newton, 
which does not fulfil the Lorentz invariance, the 
gravitational field of a spherical body does not 
depend on its state of motion, being, indeed, 
determined only by its mass M. On the con- 
trary, in the tensorial General Theory of Relativ- 
ity (GTR) by Einstein, which is a fully relativis- 
tic theory of gravitation, non-static distributions 
of matter-energy yield their own contributions to 
the overall gravitational field in addition to the 
static ones. Such peculiar terms are connected 
with the off-diagonal components goi,i = 1,2,3 
of the spacetime metric tensor, and were dubbed 
"gravitomagnetic" [1, 2, 3] in analogy with the 
magnetic fields generated by the electric cur- 
rents. Indeed, in its weak-field and slow-motion 
approximation, the fully non-linear field equa- 
tions of GTR get linearized, thus resembling 
those of the linear Maxwellian electromagnetism; 



in particular 1 , g 0i = r? i + h 0i , \h 0i \ < 1, i = 
1, 2, 3. The resulting gravitomagnetic part of the 
equations of motion of a test particle is [4] 

x % = c (hoi jk - h okji ) x k , i = 1, 2, 3, (1) 

where dots denote time derivatives, and c is the 
speed of light in vacuum. 

At great distances r from a slowly rotating 
body endowed with proper angular momentum 
S, it is 

hQi = —sAi,i = 1,2,3, (2) 
cr 

with the gravitomagnetic vector potential given 

by [3] 

A = ~^, (3) 
where G is the Newtonian constant of gravita- 



x The quantities ?7 M „, fi, v = 0,1,2,3 are the com- 
ponents of the metric tensor of the "flat" Minkowskian 
spacetime. 
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tion. The related gravitomagnetic field is [5, 1, 6] 

B g (r) = Vx A = [S-3(S-r)r]. (4) 

Among other phenomena, the resulting Lorentz- 
like, non-central acceleration due to eq. (4) 

A L t = -2 x B g (5) 

causes the secular precession of the spin a of a 
gyroscope (Schiff spin-spin effect, [7]) with fre- 
quency 

and the secular precessions of the longitude of 
the ascending node Q 

■ 2GS , , 

SZLT = c 2 a 3(l_ e 2)3/2' W 

and the argument of pericenter u 

6GS cos I . . 

WLT = - c 2 a 3 (1 _ e2) 3/2 ( § ) 

of the orbit of a test particle (Lense-Thirring 
spin-orbit effect, [8]), both in geodesic motion 
around the spinning body. The parameters 
a, e, I are the semimajor axis, the eccentricity 
and the inclination, respectively, of the particle's 
orbit. Such phenomena are often collectively 
denoted with the catchy denomination "frame- 
dragging", although also the general relativistic 
gravitoelectric de Sitter precession [9] of an or- 
biting gyroscope in the field of a static mass is a 
part of such a category [10]. 

Experimental/observational efforts have been 
dedicated in recent times to obtain empirical 
corroborations of the aforementioned predictions 
of GTR. The Gravity Probe B (GP-B) mission 
[11, 12] is an extremely refined, sophisticated and 



expensive experiment [13, 14, 15], conceived 50 
years ago, explicitly aimed to measure, among 
other things, the gravitomagnetic Schiff effect 
with four gyroscopes in a controlled environment 
enclosed in an active spacecraft orbiting the spin- 
ning Earth since 2004. The properly scientific 
phase of the mission ended in 2005, and the anal- 
ysis of the data collected during it is still on- 
going [16, 17, 18]. The expected accuracy was 
originally 1% or better, but it seems that the oc- 
currence of some unexpected systematic errors 
[15, 19, 20] may finally undermine the actual at- 
tainment of such a goal. At present, according 
to the official mission's website 2 , the claimed sta- 
tistical error is 14%, while the systematic uncer- 
tainty is 10%. 

Attempts to measure the Lense-Thirring ef- 
fect were proposed, and in some cases imple- 
mented, with some non-dedicated artificial and 
natural satellites in the Solar System; for a re- 
cent, comprehensive overview see, e.g., Ref. [21] 
and references therein. Concerning the per- 
formed analyses, the first tests date back to the 
mid of 90s [22, 23, 24, 25]; they were conducted 
in the gravitational field of the Earth with the 
non-dedicated LAGEOS and LAGEOS II geody- 
namic satellites 3 continuously tracked with the 
Satellite Laser Ranging (SLR) technique [26] by 
looking at the nodes of both the satellites and 
the perigee of LAGEOS II. Such attempts are 
still ongoing [27, 28, 29, 30, 31, 32] by retaining 
only the nodes of LAGEOS and LAGEOS II. 



2 See http://einstein.stanford.edu/ on the WEB. 

3 They are dense, spherical targets entirely covered 
with retroreflectors for passively bouncing back the laser 
impulses sent to them by ground-based stations. Both the 
LAGEOS spacecraft orbit at altitudes of about 6000 km, 
so that they do not suffer macroscopic orbit decay due to 
the atmospheric drag. As a consequence, their lifetime is 
evaluated to be of the order of 10 5 yr. 
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The claimed accuracy in such more recent tests 
is 10 - 15% [27, 28, 29, 30, 31, 32], but other 
evaluations, dealing with certain sources of sys- 
tematic errors in a more conservative way, point 
towards figures which may be up to 2 — 3 times 
larger: for such critical views, see Refs. [33, 21] 
and references therein. 

In this paper we want to clearly point out 
some epistemological and physical issues per- 
taining the performed LAGEOS-based analyses 
which have not yet been explicitly addressed in 
a satisfactorily way. In Section 2, after a brief 
review of the status of the GP-B mission and the 
perspectives for performing other measurements 
of the Schiff precession with artificial and natu- 
ral probes, we discuss if independent tests of the 
Lense-Thirring effect really exist in literature af- 
ter about 15 years since the first attempts were 
implemented. Section 3 is devoted to the rela- 
tion among the general relativistic gravitomag- 
netic field of the Earth and the spacecraft-based 
global solutions for the classical part of the ter- 
restrial gravitational field produced so far. Some 
alternative approaches to process the data of the 
LAGEOS and LAGEOS II spacecraft are dis- 
cussed in Section 4. The issue of the actual level 
of cancelation of the corrupting bias due to the 
classical quadrupole mass moment of the rotat- 
ing Earth in the tests performed so far is tackled 
in Section 5; in it the impact of the other multi- 
poles of the terrestrial gravitational field accord- 
ing to the first models from GOCE is discussed 
as well. Section 6 summarizes our findings. 

2 Do really independent tests 
of frame-dragging exist? 

Physics is an activity whose results are con- 
sidered as (provisionally) established if the ex- 



periments/observations which yielded them have 
been subsequently repeated by different teams of 
independent researchers in different laboratories 
with different methodologies. Actually, this is 
not (yet ?) the case for gravitomagnetism. 

2.1 The Schiff effect and the GP-B ex- 
periment 

Concerning GP-B, if, on the one hand, the anal- 
ysis of the data collected in 2004 — 2005 could, 
both in principle and in practice, be repeated 
by other independent researchers, on the other 
hand it will likely not be possible to do that for 
the entire experiment in any foreseeable future 
in view of its extreme sophistication and cost. 

This is certainly not satisfactorily from an 
epistemological point of view because GP-B 
seems destined to remain a unique empirical 
check of the Schiff prediction. Indeed, a proposal 
to use spacecraft orbiting the Sun and Jupiter 
[34] had not sequel so far; on the other hand, 
its complexity, cost and technological difficulties 
would certainly not have been lower than those 
of GP-B itself. 

Moreover, independent measurements of the 
Schiff effect with natural bodies in, e.g., the So- 
lar System are in all probability unfeasible. To 
this aim, let us recall that the maximum value 
of the gravitomagnetic Schiff precession occurs 
when the angular momentum S of the central 
source and the precessing spin er of the gyroscope 
are mutually perpendicular, being, instead, zero 
when they are aligned [35]. In principle, a nat- 
ural scenario satisfying such a requirement is 
the Sun-Uranus system. Indeed, while the solar 
equator is inclined to the mean ecliptic at the 
epoch J2000 by the Carrington angle Iq = 7.15 
deg [36], the spin er^ of Uranus is tilted to the 
ecliptic by 97.77 deg [37]. Of course, apart from 
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the difficulties of devising some effective meth- 
ods for continuously monitoring the precessional 
motion of the spin of Uranus, the magnitude of 
the Sun-Uranus Schiff precession would be in- 
significantly small. In principle, another poten- 
tial natural laboratory may, perhaps, be the dou- 
ble pulsar PSR J07373039A/B [38, 39]. Indeed, 
while the spin 4 Sa of A is perpendicular to the 
orbital plane [42], <tb is not aligned with Sa 
[40, 41] because of the de Sitter precession [9] 
which has recently been measured with a 13% 
accuracy [40, 41]. Actually, the gravitomagnetic 
Schiff-like spin precession [35] of <tb caused by 
Sa would be much smaller and quite difficult to 
measure. 

2.2 The LAGEOS-based tests 

The situation is, in principle, more favorable for 
the Lense-Thirring tests performed with the LA- 
GEOS satellites. 

The first attempts to reveal the existence of 
the Earth's gravitomagnetic field by analyzing 
the data of LAGEOS and LAGEOS II are 14-15 
years old, dating back to 1996-1997 [22, 23, 24]. 

The network of the 5 International Laser Rang- 
ing Service (ILRS) [43] consists of a large set 
of laser ranging stations disseminated through- 
out the world, so that the SLR community is 
quite numerous. The LAGEOS satellites, which 
are some of the most important SLR targets, 
are continuously tracked since long time. The 
GEODYN software [44], developed by NASA, 
is widely disseminated throughout the SLR sta- 
tions also because it is free of charge. More- 

4 Since the rotational periods of A and B are 23 ms 
and 2.8 s, respectively, Sa is larger than ctb- The latter 
one describe a full precessional cycle in 75 yr because of 
the general relativistic de Sitter precession [40, 41]. 

5 See http://ilrs.gsfc.nasa.gov/ on the WEB. 



over, some institutions developed their own or- 
bit analysis systems like UTOPIA by the Center 
for Space Research (CSR) of the University of 
Texas, and the Earth Parameter and Orbit Sys- 
tem (EPOS) by GeoForschugsZentrum (GFZ). 

Despite this situation, potentially favorable 
for performing several truly independent tests 
of such a prediction of GTR, none has been 
actually either performed so far, or published 
in international peer-reviewed journals. Indeed, 
apart from a pair of conference talks given by 
J. Ries et al. (CSR) [45, 46], all the relatively 
more accurate 6 tests published so far in peer- 
reviewed papers or edited books have I. Ciufolini 
in the authorship as first author or editor himself 
[27, 28, 29, 30, 31, 32]. Moreover, to date, no in- 
dependent works on the Lense-Thirring effect at- 
tributable to members of GFZ exist in literature. 
Thus, although the list of co-authors of the pa- 
pers by Ciufolini has often changed and in some 
of the most recent works [31, 32] unpublished re- 
sults obtained with UTOPIA and EPOS are de- 
scribed as well, such tests cannot be considered 
as truly independent ones. This is particularly 
true also in view of the fact that the method- 
ology adopted is basically the same, apart from 
the orbital processors used. This point will be 
explained better in Section 3 and Section 4. 

Thus, it will be possible to speak about gen- 
uinely independent tests of the Lense-Thirring 
effect with the LAGEOS satellites only if and 
when papers published in peer-reviewed jour- 
nals without Ciufolini in the authorship, and au- 
thored if possible by different researchers with re- 
spect to those having more or less systematically 
co-operated with him, will appear in literature. 



One of the major critical points of the earlier tests 
was the use of the perigee of LAGEOS II [47], heavily 
perturbed by several non-gravitational classical forces. 
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Moreover, and, perhaps, most importantly, also 
the methodology used should be different from 
that adopted so far. Such studies, which should 
be made publicly available, would be of great 
importance even in the case of negative and/or 
inconclusive outcomes. Finally, somebody may 
likely wonder why the author of the present pa- 
per does not undertake himself the task that he is 
suggesting to others. It may be pointed out that 
if, on the one hand, the results presented by skil- 
ful and experienced researchers in satellite data 
processing have raised doubts until now, on the 
other hand it is likely that analogous uncertain- 
ties would be even stronger in the case of a work 
produced by a scientist not yet actively engaged 
in such a difficult art. Moreover, the reliability of 
such results may, perhaps, be reduced in the eyes 
of a part of the community in view of the fact 
that the author of the present paper would dif- 
ficultly be considered as sufficiently neutral and 
detached from the subject considered. Analy- 
ses by really independent third parties may have 
more chances to be accepted without some sorts 
of prejudices. 

Concerning the non-negligible role played by 
such considerations, it maybe instructive to illus- 
trate the following case. In late 2007 a preprint 
titled "A critical analysis of the GP-B mission. I: 
on the impossibility of a reliable measurement of 
the gravitomagnetic precession of the GP-B gy- 
roscopes", authored by G. Forst, was posted on 
the arXiv repository [48]. This author never ei- 
ther posted other preprints on the arXiv website 
or published any peer-reviewed papers. More- 
over, there is no mention at all on the WEB of 
the organization quoted as his affiliation. Fi- 
nally, the references cited by G. Forst did not 
actually show what was attributed to them in 



the main text of his preprint, as noted first by 
K. Krogh. In early 2008, the arXiv moderators 
removed the preprint by G. Forst with the fol- 
lowing motivation: "This submission has been 
removed because 'G. Forst' is a pseudonym of a 
physicist based in Italy who is unwilling to sub- 
mit articles under his own name. This is in ex- 
plicit violation of arXiv policies. Roughly sim- 
ilar content, contrasting the relative merits of 
the LAGEOS and GP-B measurements of the 
frame-dragging effect, can be found in pp. 43-45 
of [30]." Even so, in late 2008 Ref. [48] was cited 
by I. Ciufolini in some talks of him [49, 50, 51]. 

At a different level of relevance, we also men- 
tion the 8 editing-war which involved the voice 
"frame-dragging" on Wikipedia in 2006-2007. It 
mainly consisted of the systematic and reiterated 
removal by Italian IPs followed by their almost 
immediate reinstatement, of all and only the ref- 
erences by the author of the present paper on 
some critical aspects of the Lense-Thirring tests 
with the LAGEOS satellites. On the contrary, 
the references by I. Ciufolini were never removed. 

3 New global Earth's gravity 
field solutions 

A distinctive feature of all the global Earth's 
gravity field solutions [52] obtained so far 
by several independent institutions from the 
data of the dedicated spacecraft CHAMP 
[53], GRACE [54] and GOCE [55] is that 
GTR was never explicitly solved for along 
with, say, the even zonal harmonic coefficients 
C^,0)^ = 2,4,6,... of the geopotential. The 

7 See on the WEB http:/ /www. physicsforums.com/show- 
thread.php?t=104694&page=18, post #282. 

8 See http://en.wikipedia.org/wiki/Talk:Frame- 
dragging#Recent_controversy on the WEB. 
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first global solution, obtained from CHAMP, 
dates back to 2001 [56]. Both the previous ones 
and all those of the CHAMP / GRACE /GOCE 
era are publicly available on the Internet 
at the official website of the International 
Centre for Global Earth Models (ICGEM), 
http://icgem.gfz-potsdam.de/ICGEM/. This 
fact implies, among other things, that the even 
zonals may retain a sort of a-priori "imprint- 
ing" of the Lense-Thirring effect itself; sim- 
ilar arguments were put forth in the pre- 
CHAMP/GRACE/GOCE era in Refs. [47, 57]. 
In Ref. [58] it was explicitly shown that this 
may actually be just the case for GRACE, given 
the present-day level of accuracy in estimating 
the even zonal coefficients and the size of the 
gravitomagnetic effect on the orbit of GRACE. 
Although likely time consuming, producing new 
global Earth's gravity field solutions by explicitly 
solving for relativity as well would be a really im- 
portant and independent test of the general rela- 
tivistic gravitomagnetic component of the field of 
the Earth, also because it would allow to inspect 
the correlations among the estimated solve-for 
relativity parameter(s) and the even zonals in 
the covariance matrix [47]. It would be impor- 
tant to judge if it made sense to employ that 
particular gravity model in processing the data 
of the LAGEOS satellites to try to safely ex- 
tract the Lense-Thirring effect. The inquiries of 
the present author to some scientists presently 
involved in the generation of the global gravity 
field solutions have, in general, received rather 
evasive answers, if any, mainly centered on the 
issue of the great computational and time ef- 
forts which would be required to re-process all 
the data sets from, say, GRACE. 



4 A different approach in pro- 
cessing the LAGEOS data 

An issue related to the previous one is the ap- 
proach followed to directly extract the Lense- 
Thirring effect from the data of the LAGEOS 
satellites, common to all the analyses performed 
so far. The directly observable quantity in 
SLR is the station-spacecraft range computed in 
terms of the two-way time-of-flight recorded by 
a ground-based clock. Actually, in all the tests 
implemented so far the gravitomagnetic effect on 
the LAGEOS ranges was never explicitly mod- 
elled in terms of one or more dedicated solve-for 
parameters to be estimated in the usually least- 
square sense [57], as done, instead, for a host 
of other parameters pertaining certain physical 
properties of the spacecraft, their orbital mo- 
tions and the Earth-fixed stations. Note that in 
the pre-CHAMP/GRACE/GOCE era the global 
Earth's gravity field solutions were produced 
just in such a way, i.e. by globally fitting long 
data records from a constellation of SLR targets, 
among which LAGEOS and LAGEOS II always 
played a dominant role, and estimating the even 
zonal harmonics as solve-for parameters. Inci- 
dentally, let us note that even in such circum- 
stances the Lense-Thirring effect was never mod- 
elled and solved- for. 

Another approach which may be followed may 
consist of not modeling the Lense-Thirring effect 
at all, and estimating in a purely phenomenolog- 
ical way some corrections AO to the node pre- 
cessions. This may typically be done for each 
orbital arcs. A similar approach was followed in 
the determination of the corrections Am to the 
standard perihelion precessions of some planets 
of the Solar System [59], which, in principle, ac- 
count for any unmodelled/mismodelled dynami- 
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cal effect like just the Lense-Thirring one. Also 
in the case of the timing of the binary systems 
hosting one or more pulsars a post-Keplerian pe- 
riastron precession lopk was phenomenologically 
estimated solve-for parameter along with 

other ones [60]. Subsequently, it was identified 
with the gravitoelectric precession predicted by 
GTR. 

Future independent LAGEOS-based tests 
should try to implement such strategies. In do- 
ing that, those global solutions which will be in- 
tentionally produced without retaining a-priori 
"imprinting" of relativity itself should be used 
as reference gravity field models (see Section 3). 

5 Issues pertaining the bias 
due to the Newtonian multi- 
poles of the terrestrial gravi- 
tational field 

5.1 The cancelation of the first even 
zonal harmonic of the geopoten- 
tial 

All the most recent tests performed so far rely 
upon a linear combination of the nodes of LA- 
GEOS and LAGEOS II purposely designed to 
cancel out the impact of the first even zonal har- 
monic coefficient J2 = —\/E C20 of the multi- 
polar expansion of the Newtonian gravitational 
potential of the Earth, which is a major source 
of systematic uncertainty. Indeed, it turns out 
that the nominal values of the competing secu- 
lar node precessions of LAGEOS and LAGEOS 
II caused by J2 are 7 orders of magnitude larger 
than the gravitomagnetic ones. Such a combina- 
tion was explicitly worked out in Ref. [61], fol- 
lowing a strategy put forth in Ref. [62]; see also 
Refs. [47, 63, 64]. 



Actually, the coefficient c\ of such a combi- 
nation is a function of the semimajor axes a, 
the eccentricities e and the inclinations / to the 
Earth's equator of characterizing the orbits of 
both the LAGEOS satellites. This implies that 
the unavoidable uncertainties in the computa- 
tion of such Keplerian orbital elements from the 
estimated state vectors of the satellites yield an 
overall uncertainty in c\ itself which, thus, can 
be known only with a certain numbers of signif- 
icant digits [65]. In turn, this fact introduces a 
further source of systematic bias because, given 
a certain uncertainty 8c\ w 1 x — 8 x 10~ 9 
depending on the level of accuracy with which 
one assumes in determining the inclinations I of 
LAGEOS and LAGEOS II, the resulting cance- 
lation of the J 2— induced node precessions is nec- 
essarily not perfect, contrary to what implicitly 
assumed so far. It turns out that the residual J2 
signature would amount to 14—23% of the Lense- 
Thirring one [65]. It should be remarked that a 
value of ci known up to the 8 — 9th decimal digits 
should be used to obtain just the aforementioned 
level of accuracy in the cancelation of the effect 
of J2. Instead, c\ has always been treated so far 
with a very limited number of decimal digits; in, 
e.g., Ref. [31] they are just 3 (c\ = 0.545). 

5.2 The bias due to the other even 
zonal harmonics of higher degree 

Concerning the other even zonals of higher de- 
gree Ji = —\/2£ + 1 Cefi,£ = 4,6,8, which are 
not canceled out by the aforementioned linear 
combination, their mismodeling induces a sys- 
tematic uncertainty in the resulting signal which 
may amount to a non-negligible fraction of the 
Lense-Thirring signal. The realistic evaluation 
of such a systematic alias was evaluated in sev- 
eral papers; see, e.g., Refs. [33, 21] and references 
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therein. 

Let us briefly recall that nowadays we have at 
our disposal several estimated values of the even 
zonals of the geopotential produced by different 
institutions with variable approaches and tech- 
niques from the data collected by the CHAMP, 
GRACE and GOCE dedicated missions. In eval- 
uating the aliasing of the even zonals on the 
Lense-Thirring signal it should be clear that, 
since we are dealing with the same physical 
quantities simply measured with different tech- 
niques, there are a-priori no reasons, in princi- 
ple, to prefer just one specific solution instead 
of other ones, unless objective and quantitative 
arguments are provided for trusting just it. Cer- 
tainly, in the framework of a test of fundamental 
physics, it is not acceptable to pick-up just this 
or that particular Earth'gravity models that, for 
some reasons, yield the best result in terms of fit- 
ted straight line 9 , and evaluating the systematic 
error on the Lense-Thirring effect by only using 
such particular solutions; it would be a sort of 
selection bias towards that solution just yield- 
ing the closest outcome to the one expected in 
advance. Instead, it is much more realistic and 
conservative to take into account a large num- 
ber of gravity models, provided that they are 
roughly of comparable accuracy, and adopt the 
differences among their estimated values for each 
even zonals as representative of the realistic un- 
certainty in them. In any case, quantitative, sta- 
tistical arguments should be used to discard one 
or more determinations of a given even zonal, as 
discussed in Ref. [21]. Stated simply, it is not 
admissible to play with the various gravity mod- 
els by retaining only those convenient to the a- 
priori, desired outcome and discarding, instead, 



9 The issues previously discussed in Section 3 and Sec- 
tion 4 should, at this point, not be forgotten. 



those yielding less favorable results. 

In this respect, here we point out that the first, 
preliminary global solutions from GOCE like 10 
GOCO01S, which combines data from GRACE 
and GOCE, still present significative discrepan- 
cies with respect to earlier GRACE/CHAMP- 
only global solutions. Indeed, it can be shown 
that the pair 11 GOCO01S-EIGEN51C, where 
EIGEN51C [66] is a global solution consisting 
of 6 years of CHAMP and GRACE data and the 
DNSC08 global gravity anomaly data set, yields 
an uncertainty of 27% of the Lense-Thirring sig- 
nal. Another similar example is given by the pair 
GOCO01S-AIUB-GRACE02S yielding an uncer- 
tainty as large as 23% of the Lense-Thirring 
signature; AIUB-GRACE02S [67] is a tide-free 
GRACE-only based solution obtained from al- 
most 2 yr of GRACE data. 

6 Summary and conclusions 

Whatever the final outcome of its data analy- 
sis will be, the Gravity Probe B mission, explic- 
itly dedicated to measure the general relativistic 
gravitomagnetic Schiff spin-spin effect in an ex- 
tremely sophisticated and expensive controlled 
experiment carried onboard a spacecraft orbiting 
the Earth, will remain the only empirical check 
of this specific prediction of the General Theory 
of Relativity because of the practical impossi- 
bility of repeating it in any foreseeable future. 
Moreover, no other natural laboratories in as- 
tronomical scenarios can likely be used to put 
on the test the Schiff gyroscope precession. 

See http:/ /portal, tugraz. at/portal/page/por- 

tal/TU_Graz/Einrichtungen/Institute /Homepa- 
ges/i5080/forschung/GOCO/ on the Internet. 

11 They both use the tide- free system and the fully nor- 
malized norm, so that it makes sense to compare them. 
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In principle, the situation for the tests of the 
Lense-Thirring spin-orbit effect performed so far 
in the gravitational field of the Earth with the 
non-dedicated LAGEOS and LAGEOS II satel- 
lites tracked with the Satellite Laser Ranging 
technique is more favorable. Indeed, the life- 
time of such orbiting laser targets is of the order 
of 10 5 yr, they are totally passive not requir- 
ing active instrumentation carried onboard, and 
the laser ranging community is made of several 
teams disseminated in a wide network of ground 
stations mainly using an orbital processor sys- 
tem which is freely available. Instead, despite 
the first attempts were made about 15 years 
ago by a group led by I. Ciufolini, no really in- 
dependent tests have been published so far in 
peer-reviewed journals by authors different from 
the aforementioned Italian scientist, apart from 
a couple of conference talks by a group led by 
J. Ries. On the contrary, fake Internet-based at- 
tempts to undermine the credibility of the Grav- 
ity Probe B mission were undertaken by an Ital- 
ian physicist. 

New Earth's global gravity field solutions in 
which the General Theory of Relativity is ex- 
plicitly solved for along with the multipolar co- 
efficients of the classical part of the geopoten- 
tial should be produced. On the contrary, all 
the models obtained so far from the dedicated 
CHAMP, GRACE and GOCE missions by sev- 
eral independent institutions since 2001 may be 
a-priori "imprinted" by the General Theory of 
Relativity itself since no relativistic effects were 
ever explicitly estimated in them. 

A closer connection between the gravitomag- 
netic effects on the orbits of LAGEOS and LA- 
GEOS II and the directly observable quanti- 
ties in Satellite Laser Ranging should be eluci- 
dated. In this regard, the Lense-Thirring effect 
should be explicitly modeled in the dynamical 



force models of the LAGEOS satellites, and a 
dedicated parameter should be solved-for, as it 
is common practice in all other areas of space 
science and gravitational physics. 

The cancelation of the first even zonal har- 
monic coefficient of the classical multipolar ex- 
pansion of the terrestrial gravitational potential, 
which is of degree 2, from the linear combination 
of the nodes of LAGEOS and LAGEOS II used 
so far is not perfect because of the uncertainties 
in their orbital parameters. It turns out that 
the uncanceled effect of the Earth's centrifugal 
oblateness is as large as 14 — 23% of the Lense- 
Thirring combined signal. 

Significative discrepancies among the esti- 
mated values of the even zonal harmonics in the 
first, preliminary models from GOCE and the 
earlier models from CHAMP and GRACE exist; 
according to them, the systematic uncertainty 
caused by the mismodeling in the even zonals of 
degree higher than 2 on the Lense-Thirring sig- 
nature of LAGEOS and LAGEOS II is still as 
large as about 20%. 
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